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ABSTRACT 

The recent distance determination allowed precise estimation of the orbital parameters of Cyg X-1, 
which contains a massive 14.8 M© BH with a 19.2 M© O star companion. This system appears to be 
the clearest example of a potential progenitor of a BH-NS system. We follow the future evolution of Cyg 
X-1, and show that it will soon encounter a Roche lobe overflow episode, followed shortly by a Type Ib/c 
supernova and the formation of a NS. It is demonstrated that in majority of cases (^ 70%) the supernova 
and associated natal kick disrupts the binary due to the fact that the orbit expanded significantly in 
the Roche lobe overflow episode. In the reminder of cases (^ 30%) the newly formed BH-NS system is 
too wide to coalesce in the Hubble time. Only sporadically (~ 1%) a Cyg X-1 like binary may form a 
coalescing BH-NS system given a favorable direction and magnitude of the natal kick. If Cyg X-1 like 
channel (comparable mass BH-0 star bright X-ray binary) is the only or dominant way to form BH-NS 
binaries in the Galaxy we can estimate the empirical BH-NS merger rate in the Galaxy at the level of 
~ 0.001 Myr~^. This rate is so low that the detection of BH-NS systems in gravitational radiation is 
highly unlikely, generating Advanced LIGO/VIRGO detection rates at the level of only ^ 1 per century. 
If BH-NS inspirals are in fact detected, it will indicate that the formation of these systems proceeds via 
some alternative and yet unobserved channels. 
Subject headings: binaries: close — stars: evolution, neutron — gravitation 



L INTRODUCTION 

Estimates for the rates of gravitational radiation (GR) 
sources from coalescing degenerate binaries are typi- 
cally performed with population synthesis methods (e.g., 
Lipunov, Postnov & Prokhorov 1997; Bethe & Brown 
1998; De Bonder & Vanbeveren 1998; Bloom, Sigurdsson 
& Pols 1999; Fryer, Woosley & Hartmann 1999; Nelemans, 
Yungelson & Portegies Zwart 2001, Voss & Tauris 2003 or 
more recently Belczynski et al. 2010). These studies at- 
tempt to explain the past evolution of the given observed 
binary or stellar population and put some constraints on 
the physics of stellar/binary evolution (e.g., Valsecchi et 
al. 2010). Nevertheless, there are often critical model pa- 
rameters that are poorly constrained. 

Recently, Bulik, Belczynski & Prestwich (2011) have 
taken a different approach, examining specific binary sys- 
tems with well established parameters and investigating 
their future evolution. If a binary is chosen close to the 
end of its life (e.g., the formation of double compact ob- 
ject) such a method has potentially great predictive power 
as many unknowns relating to its prior binary and stellar 
evolution can be avoided. In particular, Bulik et al. (2011) 
considered two high mass X-ray binaries (HMXBs), ICIO 
X-1 and NGC300 X-1, and showed that these systems will 
soon form close black hole H- black hole (BH-BH) systems 
that will merge within Hubble time and produce strong 
GR signature. This provided an empirical lower limit of 
detection chances for the current GR instruments without 
direct reference to population synthesis methods. 

In this study, we consider the future evolution of one of 
the most interesting binaries known in our Galaxy: Cyg 



X-1. Recently, the distance to this system was determined 
by radio parallax and other methods (Reid et al. 2011; 
Xiang et al. 2011), allowing the basic parameters of this 
binary to be firmly established (Orosz et al. 2011). This 
HMXB hosts one of the most massive (15 M©) Galactic 
BHs in a close orbit around a massive (20 Mq) O star. 
Since the companion is in the mass range for neutron star 
formation, we have selected this system to investigate yet 
another unobserved population of potential GR sources: 
black holes + neutron star (BH-NS) systems. 

2. ESTIMATES 
2.1. The future evolution of Cyg X-1 

To evolve the system forward in time we use evolution- 
ary prescriptions incorporated in the StarTrack popula- 
tion synthesis code (Belczynski et al. 2002). The evolution 
of the system is relatively simple and we do not need any 
population synthesis tools at this point. We start off with 
the best estimate of current binary parameters; the black 
hole with the mass Mbh = 14.8 M©, the optical star with 
the mass Mopt = 19.2 M© and radius of i?opt = 16.2 R© 
and the orbital period Porb = 5.6d (Orosz et al. 2011). 

The optical companion is almost filling its Roche lobe 
and will start Roche lobe overflow (RLOF) in less than 0.2 
Myr while still on Main Sequence (see Fig. [Ij. The mass 
ratio is close to unity so we do not expect the common 
envelope evolution, but rather a stable RLOF phase (Bel- 
czynski et al. 2008b, Wellstein, Langer & Braun 2001). 
However, the mass transfer rate may reach quite high val- 
ues while the donor is moving through Hertzsprung gap 
(HG). The evolution of the system is presented in Figure^ 



Mass transfer rate is calculated using physical properties 
of the donor and the system parameters (Belczynski et 
al. 2008a). The mass accretion onto BH is calculated us- 
ing the slim disk models (e.g., Abramowicz et al. 1988; 
Ohsuga et al. 2005; Ohsuga 2007) and thus can signifi- 
cantly exceed the classical Eddington limit (sec Belczyn- 
ski et al. 2008b for details). The BH increases its mass 
to Mbh = 17.8 Mq, while the optical star loses most of 
its mass Mopt = 4.2 M© to become a massive helium core 
with a bit of H-rich envelope. Note that the majority of the 
mass lost from the donor is lost from the system (highly 
non-conservative case). The period of the system increases 
to Porb = 90d. RLOF stops as the donor decreases in size 
due to the lost of its H-rich envelope. The massive helium 
or Wolf-Rayet (WR) star, ((i?wR <; 1 R©) is well within 
its Roche lobe i?iobc ~ 60 Kq. After some wind mass loss 
(^ 0.5 Mq) and about 2 Myr after the RLOF termination, 
the WR star explodes in Type Ib/c supernova and forms 
a neutron star (NS). 

There is about 2 Mq mass loss in the supernova and 
that is not enough to disrupt the system. However, the 
pre-supernova binary is rather wide (Porb = 104d, and 
semi-major axis a = 250 Rq due to additional orbit ex- 
pansion caused by the mass loss from the WR star), and 
any large natal kick tends to disrupt the binary. In Ta- 
ble [1] we list the disruption and survival probabilities for 
two assumptions about the distribution of kick velocities. 
The "full" kicks are adopted from the velocity distribu- 
tion of the Galactic single pulsars that is best described 
by a Maxwellian with ID cr = 265 km s^^ (Hobbs et al. 
2005) and uniform distribution of orientation. As there 
may be some observational and theoretical evidence that 
natal kicks are smaller in close binary systems (see the dis- 
cussion in Belczynski et al. 2010) we also explore a "half" 
kick model with kicks drawn from the same distribution 
but with a = 132.5 km s~^. As expected for such a wide 
system, the binary disruption and formation of two single 
compact objects is most likely: 94% and 74% for full and 
half kicks, respectively. Less likely, but still quite probable 
is the formation of a wide BH-NS system: 5% and 25%. 
The least likely is the formation of the close BH-NS system 
with the coalescence time below Hubble time (13.47 Gyr): 
0.2% and 0.8%. In order for the system to form a close 
BH-NS binary, the supernova explosion must produce a 
significant kick (T4ick ~ 150 km s""'^) that sends the NS 
toward the BH. 

Additionally, we have investigated an evolutionary sce- 
nario in which the RLOF starts early in the HG. Since 
at this point the optical star has a more massive core 
as compared with the above example, it will retain more 
mass through RLOF (only the H-rich envelope is trans- 
ferred/lost). Since less mass is transferred/lost the orbital 
expansion is not as dramatic and the final orbital period 
of the binary at the supernova explosion is Porb = 61. 7d 
(a = 174 Rq). This obviously leads to higher survival 
probabilities and higher chance of the close BH-NS forma- 
tion: 0.4% and 1.4% for the full kick and half kick models, 
respectively. 

2.2. Rate Estimates 

Cyg X-1 has been detected because of its unusually 
strong X-ray brightness. Let us consider the question of 
the formation rate of Cyg X-1 like binaries in our Galaxy. 



The upper limit on the X-ray active phase in Cyg X-1 is 
set by the evolutionary time of the secondary, which is 10 
Myr for a 20 Mq star. We employ this upper limit. Had 
we adopted more formal estimate of 5 Myr (it takes about 
5 Myr for a massive star to form a BH that is already in 
the system and start an X-ray phase) the estimated de- 
tection rate of BH-NS inspirals in gravitational radiation 
would increase by factor of 2. Given that we see just one 
such object and assuming that the current star formation 
rate is representative across the lifetime of the Galaxy we 
can infer the formation rate of Cyg X-1 like binaries as 
one per 10 Myr, i.e. r ^ 10~^ yr^^. Our simulations 
show that the chance of forming a merging BH-NS system 
is between 2 x 10~^ and 14 x 10""^ (see bottom row of 
Tab.[T|) from a Cyg X-1 like progenitor. These means that 
the formation rate of BH-NS binaries from Cyg X-1 like 
progenitors will be only 2 — 14 x 10~^° yr~^. This means 
that only a handful (2 — 14) close BH-NS systems might 
have been formed over a 10 Gyr history of the Milky Way. 
For comparison the Galactic empirical NS-NS merger rate 
is significantly higher: NS-NS 3 - 190 x 10"^ yr^^ (Kim, 
Kalogera & Lorimer 2010). 

The implied Advanced LIGO/VIRGO detection rate fol- 
lows once we determine the range of these detector for a 
BH-NS system. Given the final mass of the BH (A/bh = 
17.8 Mq) and assuming the mass of a NS to be canonical 
1.4M Mq, we obtain the chirp mass of the newly formed 
system to be McMrp = (MBHMopt)3/^(MBH + Mopt)"^/^ = 
3.8 Mq. Since the range for the Advanced detectors for 
NS-NS system with a chirp mass of 1.2 Mq is 300 Mpc, 
we obtain the range for such BH-NS systems of 786 Mpc 

(the detection distance scales like ex M^^-). If we adopt 
the Milky Way like galaxies density in the local Universe 
to be 0.01 Mpc-3 (g g^^ O'Shaughnessy et al. 2008), then 
within such distance there should be 2.0 x 10^ galaxies. 
Combining it with the Galactic rate we obtain the detec- 
tion rate in the range 0.4 — 2.8 x 10~^ yr~^, or a detection 
every 36 to 250 years in the Advanced LIGO/VIRGO. 

The observational uncertainities on component masses 
(^ 2 Mq) do not play a significant role on our findings 
and qualitatively they do not change any of our results. 
The same applies to the stellar wind mass loss rate from 
the companion star as the 20 Mq loses only about 1 Mg 
during its main sequence (e.g., Vink, de Koter & Lamers 
2001). The major uncertainities arise in the orbital evo- 
lution during RLOF in which ^ 15 Mq is lost/exchanged 
and supernova outcome (natal kicks). 

3. DISCUSSION 

So far we have discussed only one binary, Cyg X-1, 
as a potential progenitor of BH-NS system. Tomsick & 
Muterspaugh (2010) list several known Galactic HMXBs 
with a NS and a companion massive enough to poten- 
tial produce a BH: Vela X-1 (24 Mq), XTEJ1855-026 
(25 Mq), 4U1907-h09 (28 Mq) and GX301 (40 Mq). The 
expansion of massive companions will eventually lead to 
a RLOF. These systems due to extreme mass ratio will 
evolve into common envelope that will result in a merger 
of a NS with its companion aborting the BH-NS forma- 
tion. Such outcome follows from the fact that the mass of 
a NS is so small in respect to the mass of the envelope of a 
companion that there is not enough orbital energy to eject 



companion envelope (e.g., Webbink 1984). An extragalac- 
tic system LMC X-3 hosting a 10 Mq BH was belived to 
have a companion massive enough to form a NS. However, 
recent spectral analysis that takes into account irradiation 
of the companion, indicates that it is a B5 dwarf setting 
its mass at about 5 Mq (Val-Baker, Norton & Negueruela 
2007) and makes it a WD progenitor. 

HMXBs like Cyg X-1 are wind fed, and thus their non- 
degenerate stars do not fill their Roche lobes. It is there- 
fore curious that those binaries for which good system pa- 
rameters are known tend to be close to RLOF. Cyg X-1 
is an example: our analysis implies that the system will 
begin RLOF in 0.2 Myr, after a lifetime about 50 times 
longer than that. An even more extreme case is that of 
LMC X-1, in which an 11 Mq black hole is found in a 3.9 
day orbit with a 32 Mq companion. The companion cur- 
rently fills over 90% of its Roche lobe, implying that the 
system will undergo RLOF within w Q.l Myr (Orosz et al. 
2009). Note that the high mass of the companion implies 
that RLOF will result in unstable mass transfer and the 
likely formation of a common envelope system — this LMC 
X-1 is unlikely to produce any sort of double degenerate 
binary. There are only a few HMXBs with precisely known 
system parameters, so this tendency to be alarmingly near 
RLOF may be a statistical anomaly. Alternatively, it may 
be a selection effect: a system in which the companion is 
far from RLOF will have a smaller fraction of its stellar 
wind accrete onto the compact object, and thus a lower 
X-ray luminosity. 

Nevertheless, it may be worth entertaining the idea that 
this effect is neither a statistical glitch nor an observational 
bias, but that there is some unknown physical process that 
halts the growth of the companion star near the Roche 
lobe boundary. We note that the gravitational potential 
becomes shallow near the Li point, and X-ray irradiation 



becomes a bigger effect for companion stars that present 
a relatively large cross section. Both of these effects may 
dramatically change the surface structure and stellar winds 
of the companion star as it approaches the Roche lobe, 
conceivably in a self-limiting way. Such winds might be 
observable through the strength, shape and variability of 
emission lines associated with the wind. If such an ef- 
fect is present, the evolutionary path of the binary system 
may prove to be quite different from what is commonly 
assumed, and what we have assumed above. 

To evaluate the possible effect this might have on rates 
of GR sources, we consider how the future evolution of 
an HMXB like Cyg X-1 might proceed if it never reaches 
RLOF. As a limiting case, we constrain the orbital period 
to remain at its currently observed value of Porb = 5.6d. 
In fact, we expect the orbital period to increase somewhat 
(even without RLOF) as mass is lost from the system in 
stellar wind. In this case, the survival through the su- 
pernova and the formation of the close BH-NS system is 
expected in 8.4% of cases and corresponds to Advanced 
LIGO/VIRGO detection rate of '^ 1 per decade. So de- 
spite the fact that we have violated (in favor of producing 
close BH-NS systems) our current understanding of the 
stellar evolution we still do not get enough BH-NS merg- 
ers to expect detection in gravitational waves. 

Thus we find that if indeed BH-NS mergers are observed 
as GR sources, their immediate precursors will not be sys- 
tems like Cyg X-1 that are currently observed. 
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Wave Astronomy (KB) and NSF-AST grant 0707627 
(CB). 
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Table 1 
BH-NS Formation Statistics'' 



104d 



61. 7d 



5.6d 



SN disruption 
Wide BH-NS 
Close BH-NS 



0.944 (0.743) 
0.054 (0.249) 
0.002 (0.008) 



0.914 (0.673) 
0.082 (0.313) 
0.004 (0.014) 



0.680 (0.446) 
0.237 (0.500) 
0.083 (0.054) 



** Fraction of Cyg X-1 like systems that after the second supernova wiU be 
disrupted, wiU form a wide BH-NS or close BH-NS (merger time shorter than 
the Hubble time). The fractions are given for the full natal kicks with a — 265 
km s~^ (or half kicks with a = 132.5 km s~^). 

^ Numbers for Porb = 104 and 61. 7d correspond to physical system modeling 
with RLOF starting while the optical star is on MS and HG, respectively. The 
last model is unphysical, no RLOF was assumed and the orbital period was kept 
constant at Porb = 5.6d through the evolution (see the Discussion). 
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Fig. 1. — Radius evolution of optical star in Cyg X-1. Current radius is found at P := 16.17P0 and that places the star 
at the end of its Main Sequence (dashed line) or at the beginning of Hertzsprung gap (dot-dashed line) . Since the star is 
very close to its Roche lobe (Piobc = 17.24Po for the orbital period Porb = 5.6d) and since it has not yet started RLOF it 
means that the star is on Main Sequence at the first intersection of its evolutionary track and radius line of P = 16.17P0. 
Once the star increases its radius by about IP© it will start RLOF while still on Main Sequence and the RLOF will 
last through Hertzsprung gap (see Fig. [2]). The radius evolution is taken from single star models of Hurley, Pols & Tout 
(2000). 
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Fig. 2. — Evolution of Cyg X-1 through RLOF that will start in about 10^ yr. Bottom panel: Mass transfer rate from the 
massive donor star is very high. However, it is much lower while the donor is on Main Sequence (1 — 3 x 10~^ Mq yr~^: 
dashed line) as compared to the transfer during Hertzsprung gap (10""^ — 10"'' Mq yr"-': dot-dashed line). Mass accretion 
onto BH is factor of ^ 3 — 5 lower than the transfer rate to account for the fact that BH cannot accept all the transferred 
material. Note that the accretion rate is significantly higher than typically employed Eddington rate (5 x 10~^ M© yr"-' 
for a 15 Mq non-spinning BH) as we account for more realistic super critical accretion (slim disk ADAF) onto a rapidly 
spinning BH with a = 0.9 (Gou et al. 2011). Middle panel: The donor loses most of its mass to become a 4 M© helium 
core with a small H envelope. Most of the donor mass (~ 9 Mq) is lost from the system, while BH increases its mass from 
14.8 to 17.8 Mq. Top panel: Period of the system changes from currently observed 5.6 days to 90 days. System becomes 
wide after the RLOF due to the non-conservative mass exchange and mass ratio reversal (most of the mass is accreted 
onto BH while the donor became the less massive component of the binary). 



